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Abstract. On the basis of the best available member list 
and duplicity information, we have studied the radial dis- 
tribution of 270 stars and multiple systems earlier than 
KG in the Pleiades. Five new long period spectroscopic 
binaries have been identified from the CORAVEL obser- 
vations. We have found a clear mass segregation between 
binaries and single stars, which is explained by the greater 
average mass of the multiple systems. The mass function 
of the single stars and primaries appears to be signifi- 
cantly different. While the central part of the cluster is 
spherical, the outer part is clearly elliptical, with an el- 
lipticity of 0.17. The various parameters describing the 
Pleiades are (for a distance of 125 pc): core radius rc — 0?6 
(1.4 pc), tidal radius rj = 7?4 (16 pc), half mass radius 
rm/2 = 0?88 (1.9 pc), harmonic radius r = 1?82 (4 pc). 
Low-mass stars (later than KO) probably extend further 
out and new proper motion and radial velocity surveys 
over a larger area and to fainter magnitudes would be 
very important to improve the description of the cluster 
structure and complete mass function. 

Key words: Clusters: open - Individual: Pleiades - Struc- 
ture - technique: radial velocity - binary: spectroscopic 



1. Introduction 

The unity of star cluster structures has been discussed 
about thirty years ago by Kholopov (1969) as a general- 
ization of the results obtained first on the basis of star 
counts, and later of proper motion studies, covering wide 
areas in several nearby open clusters made by Artyukhina 
in the Pleiades (Artyukhina 1969; Artyukhina & Kalin- 
ina 1970), Praesepe (Artyukhina 1966) and other nearby 
clusters. He clearly showed the existence of a core - halo 
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structure and got a tentative relation between the core 
and tidal radii. The star count methods give only statisti- 
cal results because no real knowledge of the member stars 
is obtained. Results based on proper motion studies are 
more sound, because cluster members are individually se- 
lected. 

The core - halo structure of open clusters based on 
the direct identification of individual member stars has so 
far been established for a small number of open clusters, 
among them Alpha Per (Artyukhina 1972), the Hyades 
(Oort 1979), NGC 6705 (Mathieu 1984; Solomon & Mc- 
Namara 1980), NGC 2682 (Mathieu & Latham 1986). This 
is easily explained by the difficulty of identifying the true 
members in the outer part of the clusters where the den- 
sity of cluster members is low with respect to the projected 
density of field stars, i.e. one star per square degree at 2?5 
from the center in Praesepe. Due to the considerable area 
occupied over the sky by nearby clusters with diameters 
of about 10°, the number of stars to measure for proper 
motions reaches several times 10^. However, a definitive 
understanding of the cluster structure and dynamical evo- 
lution can only be obtained through a detailed knowledge 
of the properties of the member stars. Pleiades and Prae- 
sepe are, among the nearby clusters, primary candidates 
to study their overall structure, because of the possibility 
of obtaining very detailed information on individual stars. 

Mass segregation and concentration of binary stars 
towards the cluster centers are predicted by theories 
and numerical simulations of cluster dynamical evolution 
(Spitzer & Mathieu 1980, Kroupa 1995, de la Fuente Mar- 
cos 1996), but observational evidences have so far also 
been difficult to gather because of the generally too lim- 
ited cluster surface coverage and the lack of systematic 
radial velocity survey or detection of binaries in the range 
of (y.'05 to I'.'O. This traditional gap in binary period cov- 
erage is now beginning to be filled by speckle interferome- 
try observations obtained at CHARA (Mason et al. 1993) 
and direct imaging in the near-IR with adaptive optics 
(Bouvier et al. 1997). But it is still true that radial veloc- 
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ity surveys of the upper main sequences in the Pleiades, 
Pracscpc and a Pcrsci clusters arc badly needed to im- 
prove the knowledge of binarity and orbit characteristics 
to the level reached by the surveys of the solar-type stars 
(Duquennoy et al. 1992). 

Previous studies of the structure of the Pleiades 
(Artyukhina 1969; Peikov 1990) relied on proper motions 
to detect the member stars in the outer part of the clus- 
ter, but without photometric data to enable an analysis 
of the membership in the colour-magnitude diagram, van 
Leeuwen's very comprehensive study was based on pho- 
tometric observation in the Walraven photometric system 
and independent new proper motions (van Leeuwen 1983, 
van Leeuwen et al. 1986). However, the studies of the coro- 
nas of the Pleiades (Rosvick et al. 1992a) and of Praesepe 
(Mermilliod et al. 1990) showed that accurate proper mo- 
tions, reliable photometry and radial velocities are needed 
to determine reliable membership estimates for the stars 
in the very outer part of the clusters. 

New data are now available for the Pleiades due to the 
systematic radial velocity survey undertaken by one of the 
authors (JCM) with the CORAVEL scanner on the lower 
main sequence, the speckle interferometry survey done at 
CHARA (Mason et al. 1993; Mason 1996) and the near- 
IR imaging (Bouvier et al. 1997). This unique material 
describes for the first time the binary characteristics of 
the cluster stars and offers the opportunity of studying 
the radial structure and mass distribution in the Pleiades 
with much more detail. 

The first and third papers in this series on the Pleiades 
(Rosvick et al. 1992a; Mermilliod et al. 1997) have spec- 
ified the membership of stars in the corona of this clus- 
ter. In addition to proper motions and photometry from 
the literature, we have used radial velocity to improve the 
selection. We have shown that about 50% of the cluster 
members in the spectral range F5-K0 are located outside 
the classical area studied by Hertzsprung (1947) and mem- 
ber stars have been found as far as nearly 5° from the cen- 
ter. This makes the real size of the cluster much larger. 
The improved list of cluster members bears a lot of im- 
portance in the determination of the mass function of the 
Pleiades. Searches for lower mass stars, extending those 
conducted in the central part (Stauffer et al. 1991, Jame- 
son 1993), will certainly reduce the discrepancies between 
the field star mass function and those observed in star 
clusters, usually explained by cluster evaporation. 

This paper will discuss the new mass function and the 
radial structure of the Pleiades, with emphasise on mass 
segregation for the brighter and fainter stars, or between 
the single and binary members. We first present the results 
of the radial velocity survey for the central part of the 
Pleiades and the orbit of a new binary (Sect. 2), discuss 
the binarity of the upper main sequence (Sect. 3), describe 
the catalogue of member stars (Sect. 4) and present the 
results on the cluster structure and the mass function in 
Sect. 5. 



2. Observations and Results 

2.1. CORAVEL observations 

Previous papers presented the radial velocities obtained 
for stars in the corona (Rosvick et al. 1992a; Mermilliod 
et al. 1997), but so far the data for stars in Hertzsprung 
(1947) catalogue have not been published. The observa- 
tions were obtained with the CORAVEL radial-velocity 
scanner (Baranne et al. 1979) installed on the Swiss 1- 
m telescope at the Haute-Provence Observatory (OHP) 
for stars later than spectral type F5 and brighter than 
B = 12.5 - 13.0. A few stars, mainly F5 stars, could not 
be observed because of their too large rotation. Between 
September 1978 and December 1995, five to nine observa- 
tions per star were obtained. It has not been possible to 
get recent epoch radial velocity observations for the com- 
plete sample, but a large fraction has been reobserved to 
improve the detection of longer period binaries. Indeed, 
most orbits determined so far in the Pleiades are shorter 
than 1000 days (Mermilliod et al. 1992). One more or- 
bit has been determined here for HII 1338, with a period 
close to 8 days. Five new spectroscopic binaries have been 
detected (HII 120, 263, 476, 916 and 2500) with velocity 
ranges of a few km s"^ (3.1, 4.5, 4.6, 4.1, 8.9). The ra- 
dial velocities of HII 263 and 476 show a clear drift over 
15 years of observations. Double lines have been observed 
for HII 1122 by Liu et al. (1991) on JD 2447128.7, with 
velocities of 41.8 and -79.1 km s~^ for components A and 
B respectively. Our observations show only one dip. 

The HII numbers (Hertzsprung 1947), the V magni- 
tudes, the mean radial velocities Vr, the standard errors 
e (in km s""'^), the number of measurements n, the time 
intervals AT covered by the observations, the probabil- 
ity P(x^) that the scatter is due to chance (Mermilliod & 
Mayor 1989) are collected in Table 1. The remarks SBIO 
or SB20 refer to orbits obtained by Mermilliod et al. 
(1992), PHB to binaries detected in the colour-magnitude 
diagram from the same reference, IRB to binary found 
by infra-red imaging (Bouvier et al. 1997), VB to visual 
binary. Individual data are available upon request (Jean- 
Claude. Mermilliod@obs.unige.ch). The whole CORAVEL 
dataset is currently being recalibrated for possible zero- 
point error and color effect. All the individual data will 
then be published in a comprehensive catalogue of obser- 
vations in the Pleiades. 

2.2. The orbit of HII 1338 

HII 1338 had been initially discarded because of the large 

published Vs'mi value. 110 km (Kraft 1967), which 
was an artefact due to the double-lined character of this 
star discovered by Soderblom et al. (1993). We have there- 
fore observed this star over two seasons and got a quasi- 
circular orbit (Fig. 1), well defined with 25 and 21 obser- 
vations of the primary and secondary stars respectively. 
This system proved to be interesting because the period 
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Table 1. CORAVEL data for members of the Pleiades 



HII 


Tfly 






fi 


AT P(x^) 


Remarks 


25 


9.47 


4.9 


0.9 


5 


4685 


0.440 




34 


11.99 


4.9 


0.2 


4 


3631 


0.420 




102 


10.51 


3.9 


0.2 


19 


5559 


0.052 


IRB, 3. 63 


120 


10.79 


5.7 


0.6 


5 


4765 


0.000 


SBl, new 


129 


11.47 


6.5 


0.2 


7 


4765 


0.053 




152 


10.73 


5.2 


0.2 


5 


4765 


0.411 




164 


9.53 


9.0 


4.2 


4 


2124 


0.000 


SBl 


173 


10.87 


4.8 


0.1 


51 


4406 


0.000 


SB20 


186 


10.51 


5.5 


0.2 


8 


3624 


0.271 


PHB 


193 


11.29 


7.6 


0.2 


7 


4765 


0.169 




233 


9.66 


5.6 


0.6 


8 


5833 


0.000 


SBl 


248 


11.02 


6.3 


0.4 


5 


4763 


0.031 




250 


10.70 


4.7 


0.5 


6 


4733 


0.000 


SB ? 


253 


10.66 


4.2 


0.8 


6 


4780 


0.039 




263 


11.54 


4.6 


0.5 


9 


5822 


0.000 


SB, new 


293 


10.79 


5.6 


0.2 


7 


5020 


0.610 




296 


11.46 


5.8 


0.4 


7 


5526 


0.000 


SB? 


298 


10.86 


4.8 


0.2 


6 


4763 


0.473 


IRB, 5. 69 


303 


10.48 


6.0 


0.5 


7 


5053 


0.001 


VB, 1'.'7 


314 


10.58 


4.3 


0.7 


10 


5020 


0.008 




320 


11.03 


5.7 


0.2 


26 


5559 


0.000 


SBIO 


338 


9.06 


9.7 


1.0 


5 


5049 


0.719 




345 


11.61 


4.7 


0.3 


8 


5557 


0.072 




405 


9.82 


3.9 


0.3 


8 


5021 


0.351 




430 


11.40 


4.6 


0.2 


6 


5526 


0.566 




470 


8.90 


6.2 


0.6 


5 


3614 


0.591 




476 


10.80 


5.4 


0.7 


7 


2958 


0.000 


SB, new 


489 


10.38 


4.9 


0.3 


9 


5107 


0.099 




514 


10.71 


4.8 


0.3 


6 


4818 


0.091 




522 


11.96 


6.3 


0.1 


22 


2925 


0.000 


SBIO 


530 


8.95 


5.6 


0.6 


6 


4376 


0.011 




571 


11.23 


5.7 


0.1 


24 


4404 


0.000 


SBIO 


627 


9.66 


5.6 


1.3 


6 


4685 


0.000 


SB ? 


659 


12.09 


5.7 


0.3 


4 


3309 


0.596 




739 


9.55 


5.8 


0.2 


12 


5561 


0.016 


PHB 


746 


11.28 


6.5 


0.3 


9 


5559 


0.000 


SB ? 


761 


10.56 


5.4 


0.1 


33 


4693 


0.000 


SBIO 


879 


12.79 


4.6 


0.4 


3 


1483 


0.034 




885 


12.06 


5.2 


0.3 


4 


4105 


0.038 


IRB, 0. 87 


916 


11.71 


5.4 


0.8 


5 


5523 


0.000 


SB, new 


923 


10.14 


6.4 


0.3 


11 


5832 


0.052 




975 


10.57 


4.9 


0.4 


26 


5831 


0.003 


PHB 


996 


10.41 


6.1 


0.3 


7 


5020 


0.156 




1015 


10.55 


5.8 


0.3 


6 


4731 


0.190 




1032 


11.22 


4.1 


0.5 


6 


5522 


0.549 




1095 


11.92 


4.8 


0.3 


4 


4105 


0.067 




1101 


10.27 


5.4 


0.3 


8 


5020 


0.463 




1117 


10.20 


7.1 


0.1 


45 


3756 


0.000 


SB20 


1122 


9.28 


5.0 


0.6 


7 


0/oU 


U.zoz 


SB2 


1124 


12.22 


5.6 


0.5 


3 


1773 


0.052 




1132 


9.42 


5.2 


1.3 


5 


5055 


0.047 




1139 


9.38 


6.4 


1.0 


5 


4685 


0.020 




1182 


10.46 


6.2 


0.3 


5 


5526 


0.956 


IRB, 1'.'14 


1200 


9.90 


5.7 


0.5 


5 


4403 


0.031 




1207 


10.47 


5.5 


0.2 


6 


4731 


0.291 




1215 


10.53 


6.0 


0.2 


7 


5751 


0.753 





Table 1. (continued) 



1111 




Vr 


c 




A7- 




Remarks 


1220 


11.74 


6.1 


0.2 


4 


4732 


0.503 




1275 


11.46 


5.4 


0.2 


6 


5526 


0.398 




1298 


12.26 


6.4 


0.4 


3 


1506 


0.079 


IRB, 1'.'18 


1332 


12.46 


5.6 


0.2 


3 


1483 


0.874 




1338 


8.69 


6.3 


0.4 


25 


773 


0.000 


SB20 


1392 


9.50 


6.7 


0.5 


4 


1505 


0.068 




1514 


10.48 


4.9 


0.3 


7 


5108 


0.168 




1593 


11.11 


6.9 


0.2 


6 


5525 


0.669 




1613 


9.87 


5.7 


0.3 


7 


4694 


0.449 




1726 


9.27 


5.8 


0.3 


8 


5462 


0.035 


VB, 0'.'57 


1766 


9.13 


6.5 


0.5 


7 


4269 


0.141 




1776 


10.92 


5.8 


0.2 


6 


5522 


0.727 




1794 


10.36 


5.7 


0.2 


7 


5108 


0.397 




1797 


10.10 


5.9 


0.4 


7 


5751 


0.066 




1856 


10.02 


6.1 


0.4 


7 


4694 


0.029 




1924 


10.34 


5.6 


0.3 


7 


5108 


0.266 




2027 


10.91 


4.4 


0.3 


40 


5881 


0.000 


SBIO, IRB, O'.'l 


2106 


11.53 


4.9 


0.2 


6 


5554 


0.899 


IRB, 0'.'32 


2126 


11.68 


5.5 


0.2 


6 


5554 


0.656 




2147 


10.85 


11.1 


1.0 


36 


5828 


0.000 


SB2 


2172 


10.44 


5.9 


0.1 


34 


4693 


0.000 


SBIO 


2278 


10.91 


4.5 


0.2 


9 


5822 


0.034 


IRB, 0'.'37 


2284 


11.35 


6.3 


1.1 


8 


5874 


0.000 


SBl 


2311 


11.35 


5.4 


0.2 


5 


5554 


0.545 




2341 


10.87 


6.7 


0.2 


5 


4728 


0.474 




2366 


11.53 


6.2 


0.2 


6 


5555 


0.299 




2406 


11.10 


6.0 


0.1 


25 


4401 


0.000 


SBIO 


2407 


12.24 


5.7 


0.1 


21 


2562 


0.000 


SBIO 


2462 


11.52 


6.4 


0.2 


6 


5555 


0.577 




2500 


10.20 


5.7 


1.4 


6 


1740 


0.000 


SB, new 


2506 


10.23 


6.1 


0.4 


7 


4694 


0.004 




2644 


11.05 


4.9 


0.2 


5 


5522 


0.153 




2665 


11.36 


7.1 


0.2 


5 


5555 


0.539 




2786 


10.31 


5.2 


0.4 


6 


5520 


0.502 




2880 


11.75 


5.4 


0.3 


6 


5555 


0.083 




2881 


11.55 


5.2 


0.2 


7 


5555 


0.187 


IRB, 0'.'08 


3096 


12.12 


5.8 


0.2 


3 


1515 


0.616 




3097 


10.97 


5.6 


0.2 


22 


5555 


0.000 


SBIO 


3179 


10.04 


6.1 


0.2 


5 


4403 


0.724 





(Table 2) is very close to the theoretical cut-off period 
for circular orbit predicted by Zahn & Bouchet (1989). 
The cross-correlation dips exhibit a rather small constrast, 
which explains the higher than usual value of the O — C 
residuals. There is a hint of a third companion, which pro- 
duces a dip visible when the primary components are at 
maximum separation in radial velocity and with a veloc- 
ity close to the cluster value. Further observations will be 
done to clarify this point. 
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Table 2. Orbital elements of HII 1338 



Element 


Values 




Errors 


p[d] 


7.7570 


± 


.0005 


1 [HJD-2440000J 


8999.5 


± 


.3 


e 


0.035 


1 

± 


.009 


n — ii 

7 [km s J 


6.3 


± 


.4 


CJ [°] 


297. 


± 


13. 


K i [km s J 


59. 1 


± 


.0 


K2 [km s ] 
mi sin^i [Mq] 


65.2 


± 


.8 


0.82 


± 


.02 


m2 sin^i [Mq] 


0.75 


± 


.02 


ai sin i [Gm] 


6.36 


± 


.07 


a-z sin i [Gm] 


6.95 


± 


.08 


a{0 - C) [km s"'] 


2.39 






Tiobs 


25 






nobs B 


21 







> 




Fig. 2. Radial-velocity observations on JD 2438302. The open 
circle represents the O — C residual for star HII 1397 from the 
orbit by Conti (1968). This point fits well the observed trend 
during this night. 
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Phase 



Fig. 1. Radial- velocity curve for HII 1338 



2.3. Mean radial velocity 

The data contained in Table 1 form unique material to 
compute the mean radial velocity of the Pleiades cluster. 
We have excluded from the list all kind of binaries without 
an orbit and HII 2027 (variable systemic velocity) and ob- 
tained a weighted mean value of 5.67 ± 0.56 (s.e.), based 
on 67 stars. We accounted for projection effects to com- 
pute these values: we used the convergent point of Rosvick 
et al. (1992a, b) to derive spatial velocities and we pro- 
jected them on the direction perpendicular to the cluster 
center. The cluster mean radial velocity is defined in the 
lAU faint {rriy > 4.3) standard system, as it has been dis- 
cused by Mayor & Maurice (1985). This standard system 
is slightly below the reference level defined by the minor 
planets, by about 0.4 km/s (Duquennoy et al. 1991). 



3. Binarity among the upper main sequence stars 

Radial velocities for the B- and A-type stars have been 
obtained by many authors, but no systematic study was 
conducted over a long time interval. Morse et al. (1991) 
showed that it is possible to reach a precision of 1 to 2 km 
s~^ for B- and A-type stars by using a cross-correlation 
technique. Liu et al. (1991), using a similar method, pub- 
lished a small number of observations (1 to 3) for most 
member stars with a spectral type earlier than F5. 

The observations obtained at 128 A/mm by Abt et al. 
(1965), who made one of the major surveys in the Pleiades, 
suflfer from systematic errors during some nights and sev- 
eral of the binaries detected have not been confirmed by 
other observers. To investigate the problem linked to Abt 
et al.'s (1965) observations, we have looked for stars that 
appear to be single in the {V, U—B) colour-magnitude dia- 
gram and constant in the radial velocity surveys by Smith 
& Struve (1944) and Pearce & Hill (1975). This defined 
a sample of 22 stars, and the nightly behaviour has been 
examined for several nights. Observation of subsamples of 
the 22 stars indicate that on some nights the average ve- 
locity of the observed stars (-10.1 to -1-15.5 km s~^) is sig- 
nificantly different from that of the Pleiades mean velocity. 
Systematic trends have been noticed, as for example, on 
JD 2438302 (Fig. 2). 

Consequently, the binary status of the early-type stars 
in the Pleiades is not so well defined and more radial- 
velocity observations would be very useful. The poor 
present state of our binarity knowledge may appear dis- 
couraging when one thinks that the first radial velocity 
observations in the Pleiades were originally obtained in 
1901-1908 with remarkable precision by Hartmann (Jung 
1914, Hartmann 1914). 
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4. The member list 



We have composed a large table collecting all 270 mem- 
bers brighter than V = 12.0 from Hertzsprung's (1947) 
catalogue and the corona samples discussed by Rosvick ct 
al. (1992a) and Mermilliod et al. (1997). Table 3 contains 
the star identification, V and B — V from BDA, the open 
cluster database (Mermilliod 1995), the X and Y rect- 
angular positions in arc minutes, the distances from the 
cluster center, the multiplicity status, remarks and the de- 
duced masses of stars and components. Table 3 is available 
only in electronic form from the Strasbourg anonymous ftp 
server (130.79.128.5). 

The individual masses of the stars have been derived by 

different techniques depending on the multiplicity status 
of the stars. For single stars, the mass has been computed 
from the B — V colour with a power law relation between 
mass and B — V derived from an isochronc calculated by 
the models of Schaller et al. (1992), for log t = 8.00 and 
z = 0.02. 



For visual binaries, the colours have been computed 

from the magnitude difference and checked by a photo- 
metric separation of the binaries in the colour-magnitude 
diagram. From the distance above the ZAMS and the lo- 
cal slope of the ZAMS, it is possible to compute reliable 
magnitudes and colours of both components. The mass 
has then been computed from the colour indices with the 
same power law as defined for the single stars. The same 
method has been used for the so-called photometric bina- 
ries, i.e. stars that are clearly members, but without any 
positive binary detection from visual observation, speckle 
interferometry or radial velocity, but still located signifi- 
cantly above the ZAMS. For single-lined binaries, the sec- 
ondary mass has been taken within the interval defined 
by the spectroscopic minimum mass (when an orbit was 
available) and the maximum mass defined by the fact that 
the star is 5 magnitudes fainter than the primary. This 
condition is adopted because a difference of 5 magnitudes 
between the primary and secondary produces an effect of 
0.01 mag on the observed magnitude of the system relative 
to the primary initial magnitude. This procedure resulted 
in a probable excess of secondary masses around 0.5 solar 
masses. The use of red colours, like V — I, would permit 
the detection of additional secondaries, because the sec- 
ondary contribution to the total flux is more important 
in the / band than in B, and improve the mass determi- 
nation. The values of the masses contained in Table 3 are 
those adopted to study the Pleiades structure and char- 
acteristics. 
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Fig. 3. Map of the Pleiades displaying all the member stars 
considered in this study. Open circles are multiple stars and 
filled circles stand for single stars. The point sizes are related 
to the star magnitudes. The solid lines represent the galactic 
longitude and latitude, and the dashed lines are the principal 
axes of the star distribution. The drawn ellipse encompasses 
95 % of the stars. North is at the top and East at the left of 
the map. 

5. Cluster structure 
5.1. Global overview 
5.1.1. Cluster flattening 

From the data collected in Table 3, we have plotted a 

chart of the Pleiades (Fig. 3) which displays all the mem- 
ber stars considered in the present study. The filled and 
open circles represent single and multiple stars respec- 
tively. The size of the circles is proportional to the star 
magnitudes. 

The difference between single and multiple star distri- 
bution is already apparent in Fig. 3. The binary stars are 
more numerous in the central part (R < 80'). Does this 
feature depend on the binary character, i.e. on the greater 
cross-section of such systems, which favours the interac- 
tions with other stars, or only on the total mass of the 
system ? This point will be addressed in section 5.2 but 
Fig. 4 offers a first answer. It represents the same chart, 
but with the point size related to the stellar masses instead 
of the magnitudes. It is evident that the star distribution 
mainly depends on the star or stellar system masses. 

Examination of density contour maps suggested that 
the cluster outer region is elliptical in form. We applied 
a multi-component analysis to the data to derive the di- 
rections, and the ratio, of the ellipse axes. We found that 
the cluster halo is flattened roughly parallel to the galac- 
tic equator, with an ellipticity e = (1 — ^) of 0.17 (Fig. 
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Fig. 4. Map of the Pleiades displaying all the considered mem- 
ber stars . The point sizes are related to the masses of the stars 
or systems. 



3). Wc; think that this effect is real and is not due to ob- 
servational bias, because all stars listed by Artyukhina & 
Kalinina (1970) as candidates from their proper motion 
survey, which extends out to 5° from the Pleiades center, 
have been included in the CORAVEL observing list. 

Such a flattening of the outer part of open clusters is 
predicted both by theory and numerical simulations (Ter- 
levich 1987). Furthermore, Wielen (1974) predicted that 
the ratios of the three orthogonal axes of the cluster, con- 
sidered as a tridimensional ellipsoid, should be 2.0:1.4:1.0. 
The larger axis is pointing towards the galactic center and 
the smaller one is perpendicular to the galactic disk. As 
the Pleiades lie at a galactic longitude of 166?6 and are 
close to us, we only observe the ratio of the second and 
third axes, namely 1.4:1.0, corresponding to an elliptic- 
ity of 0.29. However, the Pleiades are ~ 50 pc below the 
galactic plane and we observe them under an angle of 23? 5. 
We thus observe an effective axial ratio of 1.15:1.0 only, 
which corresponds to an ellipticity of 0.13. It is slightly 
less than the value we determined. It is worth noting that 
van Leeuwcn (1983) already found a marginal flattening in 
the Pleiades, along the galactic plane. Oort (1979) found 
also a flattening in the Hyades, outside a 4 pc radius in the 
expected orientation, but the value is significantly greater 
than that predicted by theory. Finally, Gieseking (1981) 
observed a flattening roughly orthogonal to the galactic 
plane for NGC 3532, which is not predicted by theory. 

Our results on the Pleiades confirm the theoretical ex- 
pectations and van Leeuwen's results. 



5.1.2. Apparent star devoided area 

Another striking feature appearing on the maps shown 
in Fig. 3 and 4 is the apparent lack of stars at the west- 
ern borderline of the Hertzsprung area, at X ^ -50'. This 
peculiarity is already noticeable in the maps presented 
by van Leeuwen (1983). If we plot cluster maps for stars 
with magnitudes in the intervals V < 6, 6 < V < 8, 
8 < F < 10 and 10 < < 12.5, this feature is already 
present for 8 < V < 10 and is very well defined for the last 
magnitude interval. The first likely explanation available 
is an observational bias due to the extension of the well 
observed Hertzsprung area, which ends just before this 
western discontinuity. Nevertheless, it is surprising that 
the same effect is not so clearly observed on the three other 
sides of the Hertzsprung area. A more speculative expla- 
nation could be suggested by the work of White & Bally 
(1993). Following their conclusions, the Pleiades are un- 
dergoing a supersonic encounter with an interstellar cloud, 
approaching the cluster from the west. Thus it could be 
possible that the west part of the cluster suffers from a 
greater amount of extinction and therefore that this area 
appears less populated. However, the inspection of photo- 
metric data from the proper motion survey of Schilbach 
et al. (1995) does not allow us to confirm this prediction: 
the colours of the faint stars (down to magnitude V = 
18) located within the star devoided area are not different 
from those in other parts of the cluster. 



5.2. Mass segregation 

The concentration of multiple stars relative to single ones, 
of bright stars relative to fainter ones, and of massive stars 
relative to less massive ones is apparent in Figs 3 and 4. 
The segregation is even more evident in Figs. 5 and 6, 
which show a very clear increase of the star minimum 
magnitudes, or a decrease of the maximum masses respec- 
tively, with increasing radius. The stars or systems more 
massive than 2.5 Mq are all, except one, contained within 
a radius of 60', while stars or systems with masses smaller 
by a factor of two are within a radius of 180', i.e. a fac- 
tor 3. These diagrams provide more information on the 
radial distribution of stars within a cluster than the usual 
radial distribution represented in the form of a histogram 
of the number stars at various distances from the cluster 
center. In particular they show a nice relation between the 
star luminosities or masses and their distances from the 
center which could help to estimate the cluster radius at 
a given stellar luminosity or mass. Figures 5 and 6 also 
characterize the completeness of our sample in terms of 
magnitude, mass and radial extension. A similar figure, 
V magnitude versus radial distance, made by combining 
the various recent surveys which extends to much fainter 
magnitudes (Stauffer 1996, Fig. 1) shows the distribution 
of members and candidates in the Pleiades. Large gaps are 
apparent and our knowledge of the Pleiades remains in- 
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Fig. 5. Apparent magnitudes of the stars as a function of 
their radial distances to the cluster center. The open and filled 
squares denote multiple and single stars respectively. 
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Fig. 6. Logarithm of the star masses as a function of their 
radial distances to the cluster center. Same symbols as in Fig. 
5. 



complete below KO in terms of both magnitudes and area 
surveyed. 

To investigate more accurately the radial distribution 
of the different star populations we have spht the sample 
into two classes, according to different criteria: magnitude, 
mass and multiplicity of the stars. 



5.2.1. Single and multiple stars 

In contrast to the mass/magnitude parameters for which 
it is difficult to find an objective value to divide the sam- 
ple, the binary parameter is more direct: either the star 
is single or it is double or multiple and we can compare, 
without ambiguity, the cumulative distributions of the two 
populations (Fig. 7). The result is quite impressive and a 
Kolmogorov-Smirnov test clearly indicates that the two 
distributions are different. However, we should keep in 
mind that the concentration of multiple stars towards the 
cluster center may be overestimated. Indeed, the detection 
of low luminosity double or multiple systems is probably 
more complete in the central area of the cluster. In the 
central region we obtained more radial velocities over a 
longer period, which increases the completeness of the bi- 
nary detection, while, in addition, near-infrared imaging 
has been conducted only in the Hertzsprung (1947) area 
by Bouvier et al. (1997). 

This would lead one to underestimate the number of 
such systems in the cluster outer part, and to slightly 
overestimate their concentration towards the cluster cen- 
ter. However, most stars resolved by Bouvier et al. (1997) 
were classified "PHB" by Mermilliod et al. (1992), i.e. 
they were identified as binaries from their position in the 
colour-magnitude diagram. Therefore, the same reasoning 
applied to the corona ensures that we missed a minimum 
number of binary stars. 

Among the multiple star population itself, a mass seg- 
regation is obvious (Fig. 6) and the distribution of multi- 
ple stars with masses greater than 2.5 Mq (23 stars) is sig- 
nificantly different from the distribution of multiple stars 
with masses between 1.4 and 2.5 Mq (52 stars). More- 
over, if we divide the multiple star population between 
"long period" binaries (visual, IR imaging, speckle and 
occultation binaries; 44 systems) and "short period" bi- 
naries (spectroscopic binaries, 43 systems) wc notice that 
their radial distributions are very similar. These results 
are consistent with the hypothesis that the radial segre- 
gation towards the cluster center depends mainly on the 
mass of the systems, and not on their periods. The same 
conclusion was obtained by Raboud & MermiUiod (1994), 
and some theoretical considerations may help to under- 
stand these results. 

Mathieu (1985) estimated that the period boundary 
between "hard" binaries, which are systems with a bind- 
ing energy much greater than the mean kinetic energy of 
the stars in the cluster, and "soft" binaries lies at 10^-10'' 
yr in the Pleiades, if the two binary components have a 
mass of 1 Mq . Such a period implies an apparent separa- 
tion of 1'.' - 4'.'7 and, as we have only 9 stars displaying this 
characteristic, we could consider that our sample is mainly 
composed of hard binaries (78 systems among 87). Fol- 
lowing Hills & FuUerton (1980) the hard binaries increase 
their binding energies, during close encounters with sin- 
gle stars, at the same average rate independently of their 
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Fig. 7. Cumulative distributions for the multiple stars (open 
squares) and the single stars (filled squares). 

semi-major axes. Thus one cannot hope to detect any con- 
centration among the binaries in function of the period. 
But we will observe a mass segregation between the bina- 
ries, resulting from energy equipartition, because the mass 
of the multiple systems are different. However, Spitzer & 
Mathieu (1980) show that during close encounters between 
two hard binaries the softer is disrupted. This scenario is 
only valid if one of the binaries is much harder than the 
other. In that case wc should observe a greater proportion 
of hard binaries in the center of the cluster, where close 
binary encounters are most probable. Abt (1980) shows 
some evidence of such an effect, but we cannot confirm 
his result: our radial distributions of "long" and "short" 
period binaries cannot be considered as different. 



5.2.2. Sample subdivisions using a magnitude or a mass 
criteria 

Analysing the radial distribution of star populations de- 
fined by a magnitude or mass criteria implies the selection 
of a cut-off value. We could define this value as the one 
which optimises the separation between two well repre- 
sented populations. However, such a selection is fuzzy and 
is very difficult to apply consistently to different clusters. 
Therefore, we decided to analyse the morphological modi- 
fications of the cumulative distribution functions for vari- 
ous cut-off to determine the value at which the two distri- 
butions become statistically different (between = 4.0 
and 12.0, with a step of 0.1 mag, and between masses 
of 0.73 Mq and 5.0 Mq, with a step of 0.05 Mq). Two 
distributions are considered different if the probability of 
false rejection of the null hypothesis (that the two distri- 



butions are the same) is less than 5%, according to the 

test of Kolmogorov-Smirnov. 

Figure 8 represents seven of the 81 sample subdivi- 
sions considered, for the magnitude criteria. These seven 
diagrams are sufficient to reveal the expected continuous 
behaviour of the two cumulative distributions for the var- 
ious cut-off values. Interestingly, the two cumulative func- 
tions are considered as significantly different already for 
a cut-off value as great as m„ = 10.9, corresponding to 
masses M < 1 Mq. 

The same procedure has been applied to the sample di- 
vided according to the mass criterion. We found that the 
two cumulative distributions are significantly different for 
all cut-off values M > 0.8 Mq. This limit is smaller than 
the previous one. This difference results from the fact that, 
due to the high percentage of multiple stars in the central 
part of the cluster, the division of the sample in terms 
of magnitude underestimates the real mass difference be- 
tween the two populations. 

The adopted procedure reveals that the definition of 
a universal cut-off value between two populations based 
on the statistical signifiance of the difference between the 
distributions is hardly applicable, because the populations 
appear to be different even for a mass limit as low as 0.8 

To characterize the degree of mass segregation among 
different populations in the cluster, it was better to subdi- 
vise the sample in more than two groups. Fig. 9 represents 
the cumulative distributions for 4 different mass intervals. 
The Kolmogorov-Smirnov test indicates that the three 
more massive populations are significantly different, but 
that the distribution of the M < 1 Mq and 1 < M < 1.6 
Mq populations could not be considered as different. Even 
the radial distribution of the stars with masses M < 0.8 
Mq could not be considered as different from the popula- 
tion of stars with 1 < M < 1.6 Mq. Consequently, mass 
segregation is less important for the lowest mass stars 
of the cluster, a result already noticed by van Leeuwen 
(1983), and also found in the Hyades (quoted by Mathieu 
1985). It could be a general phenomenon because lower- 
mass star population is the most spatially extended, due 
to energy equipartition, and therefore is the most severely 
truncated by the galactic tidal field (Mathieu 1985). 

5.3. Characteristic radii 

A way to quantify the different radial distributions is to 
compute the values of characteristic radii, such as the ra- 
dius containing half of the total number of stars (r„/2)) the 
radius containing half of the total mass of the stars {r„^/2)i 
the core radius (rg), the tidal radius (rt) or the harmonic 
radius (f). The description of an apparent star distribu- 
tion with characteristic radii supposes that the cluster has 
a spherical symmetry, which is not the case of the Pleiades, 
at least for the outer part. However, the study of a slightly 
flattened system is well achieved with the assumption of 
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Fig. 8. Cumulatives distributions for the bright stars (open squares) and the faint stars (filled squares). Cut-off values are 
indicated for the sample subdivisions. 

spherical symmetry and does not need the use of more so- 
phisticated distributions. A simple simulation shows that 
the derivation of the core and tidal radii are unaffected 
by an artificial flattening of a cluster defined by a stellar 
density law of King (1962), with an axial ratio of 1.4:1 in 
order to reproduce the theoretically predicted ellipticity. 



5.3.1. Core and tidal radii 

The core and tidal radius were obtained by fitting, through 
the observed stellar-density distribution, the empirical 
density law of King (1962), given by 

f ^ 2 
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Fig. 9. Cumulative distribution functions for different mass 
intervals: M < 1 M© (filled squares), 1 < M < 1.6 Mq (open 
squares), 1.6 < M < 2.5 Mq (crosses) and M > 2.5 Mq (filled 
triangles) . 



where /(r) is the projected number density as a function 
of radius r, fc is a constant of proportionality, Tc and rj 
are the core and tidal radii respectively. The adopted pro- 
cedure for adjusting this density law could be divided in 
two parts: the determination of the initial values of k, Vc 
and ft for the fitting process and the fitting process itself. 

For the estimation of the initial parameters we com- 
puted a great number of King's stellar density laws, for the 
cluster divided in concentric rings, with k, rc and rt chosen 
by a Monte-Carlo procedure within large "windows" . We 
then selected the set of parameters for which the function 

was minimum: 



= Y] -TT [denSobs{i) - denSmodei{i)f (2) 
1(1) 
1=1 ^ ' 



where densobsii) is the observed projected number density 
in the ring, denSmodei{i) is the King's law prediction 
of this density and cr(i) is the uncertainty from Poisson 

- , where N{i) is the star number in 



A(i) 



statistics 

the ring and A{i) is the area of this ring). 

This set of parameters was used as the initial one for 
the second part of the adjusting procedure: the fitting pro- 
cess itself. We estimated the best-fit parameters by min- 
imizing the same function, in the standard manner. 
Incorrect distributions were rejected according to the cri- 
terion defined by Lampton et al. (1976). The estimated 
parameter uncertainties correspond to 95% independent 
confidence intervals. These ranges were estimated using 
the procedure of Lampton et al. (1976). We thus obtained 
fitted values for k, Vc and rt- 



The whole procedure is applied to 16 cluster subdi- 
visions in 4 to 20 concentric rings. The final values of 
the density law parameters, reported in table 4, are the 
weighted mean of these 16 sets of values. 

The tidal radius ought to be identical for all popula- 
tions and independent of the star masses of their members. 
However, due to the difficulty to extrapolate the data to- 
wards the external parts of the cluster for the concentrated 
populations, the values of given in Table 4 are different. 
The only reliable estimations of the tidal radius are obvi- 
ously those based on the more extended populations. 

5.3.2. Harmonic radius 

The harmonic radius f is defined following the equation 
of the potential energy of the cluster O (Chandrasekhar 
1942): 



i=l j>i 



niimj 



1 GM2 

2 f 



(3) 



is the stellar 
Vij is the distance between star i and j and M is 



where G is the gravitational constant, m 
mass 

the total mass of the cluster. 

If we estimate the apparent stellar density, F{r), in 
strips of equal area covering the cluster, it is possible to 
relate the harmonic radius to this density through the re- 
lation (Schwarzschild 1954) 



_^{£Fir)dr 



CF^{r)dr 



(4) 



In our study we have computed the apparent stellar den- 
sity in 10 and 30 strips. Each density is smoothed by com- 
puting its average between 20 equidistant strips, with dif- 
ferent orientations, uniformly distributed around the clus- 
ter center. The uncertainty associated to the harmonic ra- 
dius corresponds to la. It should be multiplied by a factor 
2 to be compared with the uncertainty of Vc and rt dis- 
played in Table 4. 

5.3.3. Comparisons of the various radii 

Table 4 presents the results for the radius containing half 

of the total number of stars (?'„/2)i the radius containing 
half of the total mass of the stars the core (rc), 

the tidal (rt) and the harmonic (f) radii for different dis- 
tributions. As stated above, the subdivision criterion cho- 
sen for the bright/faint stars and the massive/less massive 
stars is not based on a statistical test. We adopt a cut-off 
value of = 9.5 because it corresponds to the limit be- 
tween the A-F and the G-K stars. At spectral types F0-F5 
occurs the transition from radiative to convective atmo- 
sphere which produces a gap (Bohm-Vitense & Canterna 
1974) best seen in colour-colour diagrams. Due to the po- 
sition of this gap in the Pleiades main sequence, it is a 
good limit to split the sample. The corresponding cut-off 
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in terms of mass is M = 1.5 Mq. It is worth noting that 

the two s\ib-samplcs created with the mass siibdivision arc 
not exactly the same as the ones created with the mag- 
nitude cut-off, because the mass subdivision takes more 
efficiently into account the multiple status of the systems 
present in the cluster. 



Table 4. Characteristic radii ['] for different populations. The 
errors associated with r„/2 and are typically between 5 
and 10 [']. 



Population 


rn/2 


rm/2 




rc 




n 




r 


Complete sample 


62 


53 


38 


(14) 


445 


(182) 


109 


(32) 


Bright stars 


49 


36 


33 


(25) 


241 


(189) 


86 


(24) 


Faint stars 


72 


69 


51 


(25) 


434 


(220) 


123 


(28) 


Massive stars 


43 


37 


22 


(14) 


327 


(273) 


75 


(21) 


Less massive stars 


77 


73 


58 


(29) 


436 


(224) 


134 


(31) 


Single stars 


81 


69 


55 


(27) 


461 


(253) 


133 


(29) 


Multiple stars 


42 


35 


32 


(20) 


176 


(79) 


70 


(21) 



As expected, the different radii in Table 4 behave sim- 
ilarly. The bright, massive and multiple populations are 
charaterised by a small radius, while the faint, less mas- 
sive and single populations are characterised by a larger 
radius. Moreover, for each considered population the half- 
mass radius is always smaller than the radius containing 
the half total number of stars. This result means that we 
observe mass segregation even in each sub-sample of clus- 
ter stars. This again indicates that the segregation among 
cluster stars is mainly governed by the mass. 

5.4- The frequency of multiple star systems 

As the multiple stars are preferentially concentrated to- 
wards the cluster center, relative to single stars, it is in- 
teresting to compare their proportions in different parts 
of the cluster with predictions from numerical models. 

Kroupa (1995) predicted the existence of a large pro- 
portion of binaries in the central region of star clusters, 
computed with the relation: 

f = m 

inaries / [^binaries ~t~ ^single stars] 

In particular, he expected that the total proportion of bi- 
naries in the central 2-pc sphere of the Pleiades is "prob- 
ably" close to 60 %. This value is larger than that we 
obtained (48 %) in the central 2-pc disk. However, as is 
usually mentioned in the discussion of binary detection, 
the observed value is most probably a lower limit. In the 
outer part of the cluster, the binary frequency is 20 % 
only. 

A binary fraction dependence on primary mass is con- 
sidered as a discriminator betweeen the "capture" binary 
formation mode and the fragmentation mode. In the first 



model the binary fraction is a strongly increasing func- 
tion of primary mass, but in the second one the binary 
fraction is only weakly dependent on primary mass, and 
then in the sense of binary fraction declining with increas- 
ing primary mass (Clarke 1997). As the detection of spec- 
troscopic binaries in the halo of the Pleiades is probably 
not complete, we restricted our study of the percentage 
of multiple systems, as a function of the primary masses, 
to the Hertzsprung area. Figure 10 represents the values 
observed among this sub-sample, indicating a possible de- 
pendence of the multiple star fraction with primary mass. 
Six among the 7 brightest cluster stars belong to multiple 
systems (Fig. 5), i.e. the most massive stars are multiple 
ones. 

Before claiming this effect has a physical origin we need 
to examine possible biases due to statistics and detection 
limits. Statistically, if we consider a star cluster only com- 
posed of binaries, which were formed by random pairing, 
we could expect a decrease of the probability to observe 
a star as a primary with decreasing mass (Kroupa et al. 
1996). In such a cluster, the most massive star has a prob- 
ability of 100 % to be observed as a primary. However, 
considering a Salpeter mass function {f{m) = Cm~^'^~^^\ 
with {l+x) = 2.35) extending down to 0.08 M©, we found 
that a 0.7 M0-star (i.e. the least massive stars in our sam- 
ple) has a probability of 95 % to be observed as a primary: 
obviously the mass range from which a secondary can be 
picked decreases with decreasing mass of the primary. Us- 
ing the mass function exponent we derived for our sample 
{{1+x) =2.5, see next section), we obtain a probability of 
96 % to observe a 0.7 M0-star as a primary. This probabil- 
ity is only 82 % if we consider the exponent values used by 
Kroupa et al. (1991) and McDonald & Clarke (1993) (2.35 
for stars with mass m > 1 M©, 1.1 for 0.08 < m < 0.5 
M0 and 2.2 for 0.5 < m < 1 Mq). 

On the other hand, the detection of low mass secon- 
daries in spectroscopic binaries is more difficult for B- and 
A-type stars. So, if any, the detection biases would lead 
to underestimate the binary frequency for star with mass 
larger than 2 Mq. 

As our frequency of multiple systems range from 
~100% for the brightest sample stars to around 40 % for 
the less massive stars, we conclude that this effect could 
not be imputed to the mentioned statistical, or detection, 
biases. This conclusion also implies that the binaries are 
not created by random pairing, if the effects of dynamical 
evolution can be neglected. 

The results presented in Fig. 10, and the fact that the 
most massive cluster stars belong to multiple systems, 
could be explained by the effect of dynamical evolution 
instead of a different formation mode. Encounters with 
exchange reactions, during which massive single stars re- 
place lower mass binary components is a likely occurring 
phenomenon (Mathieu 1985). In a 100-body equal-mass 
cluster with an 11 % population of slightly hard binaries. 
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Fig. 10. Percentage of multiple stars as a function of the pri- 
mary mass, in the cluster central part. From the left to the 
right the bins involve: 28/23 (single/multiple stars); 29/18; 
11/9; 4/6 and 3/6. The dotted lines correspond to the per- 
centage of spectroscopic binaries, relative to the total number 
of stars. 

24 crossing times are enough in order that approximately 
two thirds of the binaries had encounters (Aarseth 1975). 

The Pleiades have an age of 10^ yr and therefore do 
not present a snapshot of pure star formation products. 
Thus our results cannot be used to discriminate between 
different binary formation scenarios. We need to repeat 
the same kind of investigations in very young (and pop- 
ulous) open clusters to be able to constrain theoretical 
models. However, it is worth noting the current lack of 
such a suitable observational sample, because sufficiently 
populous young clusters are distant objects suffering from 
extreme contamination problems. 

5.5. The mass functions 

The improved material collected in this paper allows var- 
ious determinations of the Pleiades mass function, for the 
complete sample, for single stars and primaries only, and 
allows examination of the effect of radial extension. 
We fitted a Salpeter-type power law in the form 

lo9(^-^)=C-{l + x)log{M) (5) 

throughout the observed data. In Eq. (5) df/dM is the 
number of stars per unit mass as a function of mass M. C 
is a constant and (l + a;) is the power law exponent, which 
has the value of 2.35 following Salpctcr (1955). 

We first derived the mass function for the complete 
sample, without any correction for the binary content of 
the cluster. We therefore computed all the stellar masses 



with our relation between [B — V) and the mass (Sect. 
4). The mass function that wc obtained is polluted with 
unresolved binaries (case 1 in Table 5). The derived mass 
function slope (2.1 ± 0.21) agrees marginally with the 
canonical value of Salpeter (1955). 

Using the available binarity information of our sample, 
we derived a mass function only for the single stars and the 
primaries for the complete sample (case 2 in Table 5). The 
slope of 2.5 ± 0.15 is not very different from the previous 
one, indicating that the determination of mass function is 
not seriously affected by unresolved multiple stars. This 
finding confirms the results of Tarrab (1982). 

It is also very interesting to compare separately the 
mass functions of the single stars and of the primaries of 
multiple systems. Following Vanbeveren (1982), the sim- 
ilarity between these two kinds of function is rather im- 
probable. He also found, in the case where the probability 
of the formation of binaries increases with increasing cloud 
mass, that the initial mass function of the primaries is less 
steep than for the single stars. These computations were 
done with masses between 15 and 150 Mq. Although 
our mass range is very different, it is worth noting that 
our results (cases 3 and 4 in Table 5 and Fig. 11) agree 
with these conclusions. This result is expected from the 
observations presented in the previous section about the 
percentage of multiple systems as a function of the pri- 
mary masses. 

Different initial mass functions for single stars and pri- 
maries would imply that these two populations have dif- 
ferent origins, and therefore discredit the "capture" , in the 
sense of random pairing, formation mode for the binaries. 
However, as noted before, we arc not observing the initial 
conditions of stellar formation and the currently different 
mass functions of the single stars and the primaries may 
result from the dynamical evolution of the cluster. Due 
to encounters with exchange reactions, the low mass pri- 
maries would tend to be replaced by massive single stars, 
decreasing the slope of the primaries mass function rela- 
tive to the single stars mass function as it is observed. If 
we now restrict the comparison of these two mass func- 
tions to the Hertzsprung area, where the multiple stars 
detection is more complete but where the mass segrega- 
tion effect has depleted the low mass stars, we still observe 
the same results. The slope (1 + x) is 1.57 ± 0.34 for the 
single star mass function and 0.68 ± 0.16 for the primary 
mass function. 

To test the effect of an incomplete surface coverage 
of the cluster field on the power law exponent, we only 
considered the stars within the Hertzsprung area (case 5 
in Table 5). The mass function slope is then less steep than 
for the complete sample (case 2). This result is evidently 
due to mass segregation, which depletes the inner cluster 
part from the less massive stars (Pandcy ct al. 1991a, b). 
It is therefore very important to consider the cluster and 
its surrounding halo in order to derive a realistic mass 
function, as already stated by Scalo (1986). If we only 
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Fig. 11. Mass functions. The solid line stands for the complete 
sample (single stars and primaries). The long-dashed line rep- 
resents the mass function of the primaries and the short-dashed 
line stands for the single stars. 



consider the outer part of the cluster (case 6) we obtain 
a very steep slope for the mass function, also expected 
from the effect of mass segregation, as this area is highly 
depleted in massive stars. 

Finally, if we compute the mass function of the cluster 
by considering the stellar system masses merely as the 
sum of the different component masses (cases 7 and 8), 
we obtain smaller slopes for the complete cluster and for 
the Hcrtzsprung area than for the corresponding single 
and primaries mass functions (cases 2 and 5). This results 
because the mass summed method "creates" massive stars 
which have no physical reality. Nevertheless the effect of 
mass segregation is still noticeable. 



Table 5 



Values of different power law exponents (1 -|- a;) 



Sample 


(1 


+ 


x) 


(1) Complete sample (with unresolved binaries) 


2.10 


± 


0.21 


(2) Complete sample (singles -|- primaries) 


2.50 


± 


0.15 


(3) Singles 


2.75 


± 


0.22 


(4) Primaries 


1.47 


± 


0.27 


(5) Hertzsprung sample (singles + primaries) 


1.71 


± 


0.14 


(6) Cluster outer part (singles + primaries) 


3.20 


± 


0.16 


(7) Complete sample (mass summed) 


1.76 


± 


0.15 


(8) Hertzsprung sample (mass summed) 


1.22 


± 


0.18 



5. 6. Estimation of the cluster total mass 

The cluster total mass is an important parameter for the 

comparison between results from numerical simulation of 
open cluster dynamical evolutions and real objects (Wie- 
len 1975). However, its derivation is very difficult (Bruch 
& Sanders 1983) as will be shown in this section. We use 
three methods to estimate the cluster total mass: the tidal 
radius, the Virial theorem and the mass function. 

5.6.1. The tidal radius 

The total mass of an open cluster is related to the tidal 
radius through 

4^(^-B) , 



G 



(6) 



where G is the gravitational constant, rt is the tidal radius 
of the cluster, A and B are Oort's constants of galactic ro- 
tation. This equation is directly derived from King (1962), 
with the assumption that the cluster is at the same dis- 
tance from the galactic center as the Sun. 

The tidal radius considered in Eq. (6) is measured in 
the direction of the galactic center. However, we only ob- 
serve the tidal radius perpendicular to this direction but 
parallel to the galactic disk. In Sect. 5.1.1 we found that 
the Pleiades have a flattening close to the expected one. 
We could therefore consider that we have a cluster with 
axes ratios 2:1.4:1 and then the tidal radius in the direc- 
tion of the galactic center has a value of 2/1.4 times the 
value of the observed tidal radius (445' from Table 4). 

Using A = 15 km s"^, B = —12 km s""'^ and rt = 
636' = 23 pc, we derive a total cluster mass of 4000 Mq. 
The confidence interval, based only on an uncertainty of 
la for rt, is 1600 to 8000 Mq. If we do not apply any 
correction due to the cluster flattening we obtain 1400 
M0 for the mass, with a confidence interval between 530 
and 2900 M0. As may be seen from Eq. (6), this cluster 
mass determination is very sensitive to the tidal radius 
value, which is only poorly determined. 

5.6.2. The Virial theorem 

To derive the cluster mass through the Virial theorem, 
we need to estimate the velocity dispersion in the cluster. 
From 67 member stars, located in the Hertzsprung area 
(see section 2.3), we obtain a radial velocity dispersion of 
0.36 km s~^, from the quadratic difference of the total 
observed dispersion corrected for projection effects (0.56 
km s~^) and the mean measurement error (0.43 km s~^). 
Assuming velocity isotropy and taking for the harmonic 
radius r ~ 109 ['arc] = 3.96 pc, we derive a cluster total 
mass of 720 Mq, through the relation 



Mr = 



6fVl 
G ■ 



(7) 



The confidence interval, based only on the harmonic ra- 
dius uncertainty, is 510 to 940 Mq, 
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5.6.3. The mass function 

Finally, we estimate the cluster mass using the observed 
mass function (see Sect. 5.5). We assume that the mass 
function for the single and primaries (case 3 in Table 5) is 
applicable down to 0.08 M© and we include a factor 1.16 

to take into account binary companions (32 % of the stars 
are binaries with a companion mass equal to the half of 
the primary mass). The result is a cluster total mass of 
950 M0, with a confidence interval between 800 and 1150 
Mq, based only on the mass function slope uncertainty. 

5.6.4. Pleiades total mass 

The three different mass estimations (Table 6) are rela- 
tively close together, within a factor 2, except the first de- 
termination taking into account a correction for the clus- 
ter flatness, which is a factor more than 5 greater than the 
lowest mass estimation, van Leeuwen (1983) estimated a 
cluster total mass of 2000 M©, which is within our con- 
fidence intervals of mass determination through the tidal 
radius. 



Table 6. Results of the different cluster total mass determina- 
tions (see text for more details). 



Method Cluster total mass 1 a confidence interval 

Mq Mq 

Tidal radius 4000 [1600, 8000] 

(with correction for the cluster flatness) 

Tidal radius 1400 [ 530, 2900] 

(without correction for the cluster flatness) 

Virial theorem 720 [ 510, 940] 

Mass function 950 [ 800, 1150] 

Summed mciss 412 



If we sum up all the stellar masses derived for our 
whole sample of stars, we obtain 412 Mq. Then the stars 
down to ruv = 12.5 and extending out to 5° from the 
cluster center represent at best ~ 57 % and at worst ^ 10 
% of the total cluster mass. Member stars fainter than V 
= 12.5 are already known in the Pleiades and about 600 
flare stars (spectral type K2 and later) have been detected. 
Assuming a mean mass of 0.5 Mq and 400 faint members, 
not all flare stars are members (Jones 1981), we can add 
at least 200 solar masses. 

These results clearly point out the great difhculties 
of estimating the total cluster mass with precision. All 
the estimators used arc based on strong hypothesis and 
are dependent on poorly determined parameters. The only 
way to properly compute a cluster mass is to sum all the 
masses of the individual cluster members. 



6. Conclusion 

A study of the Pleiades structure has been performed on 

the basis of the presently available data which limits the 
sample to stars brighter than V = 12.5. We used the best 
present knowledge on duplicity in the Pleiades. 

Using a multi-component analysis applied to the ap- 
parent stellar positions we find an ellipticity in the cluster 
outer part, in agreement with theoretical expectations and 
van Leeuwen's results (1983). 

We have observed a clear mass segregation, which de- 
pends on the mass of the stars or systems and not on 
the binary periods. Consequently, binaries are more con- 
centrated than single stars and massive binaries are more 
concentrated than less massive ones. The mass segrega- 
tion is significant down to 1 Mq . Different radii have been 
computed to characterize the radial distributions of vari- 
ous cluster star populations. 

For the first time, to our knowledge, the mass func- 
tion of single stars and primaries of multiple systems have 
been determined separately and compared. They turned 
out to be different, in agreement with predictions made 
by Vanbeveren (1982). Such a result, if confirmed by sub- 
sequent studies, may have important implications for star 
formation models. It shows the extreme interest of de- 
tailed studies of young, or very young, open clusters and, 
especially, of their binary populations. 

We review the great difficulty in deriving an estima- 
tion of the cluster total mass. The mass estimates span 
an order of magnitude, from ~500 to 8000 Mq, once we 
consider parameter Icr errors. Considering only the stars 
included in our sample we derived a projected spatial mass 
density of 17.7 Mq pc~^, or 9.5 stars pc~^, for the 2-pc 
radius central disk of the cluster. This region was chosen 
because it is widely used by Kroupa (1995) for compar- 
isons between numerical models. Central density values 
for different clusters should be compared only if they are 
computed in the same manner, for example inside this 
central 2-pc radius disk. 

There seems to be no other alternative to determine 
the Pleiades total mass and complete luminosity function 
than identifying all members in an area even larger than 
that investigated here and to fainter magnitude. The deep- 
est surveys made in the Pleiades have been limited to the 
central region within a radius of 3°. They should be ex- 
tended to at least 6° from the center. 
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